Itineracy effects on spin correlations in ID Mott insulators 
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We consider spin correlations in the one dimensional half filled repulsive Hubbard model. For 
very large values of the on-site repulsion U the spin correlations are dominated by virtual hopping 
processes of electrons and are described in terms of a spin-1/2 Heisenberg chain. As U is decreased 
real hopping processes of electrons become important and eventually dominate the spin response. 
We discuss the evolution of the dynamical structure factor as a function of U. We comment on the 
relevance of our results for inelastic neutron scattering experiments on quasi-lD materials. 
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Quasi-one-dimensional (Q1D) Mott insulators are a 
paradigm for the importance of strong correlations and 
are known to exhibit a wide variety of unusual physical 
phenomena such as spin-charge separation and quantum 
number fractionalization. Experimental realizations in- 
clude the chain cuprates SrCuOo and Sr^CuOs and a 
number of organic compounds 0,12 ■ It has been known 
for a long time that the magnetic properties of Q1D Mott 
insulators are rather unusual. The dynamical structure 
factor is entirely incoherent and reflects the fact that 
the elementary spin excitations carry spin 1/2. They 
are very different in nature from antiferromagnetic spin 
waves. In the limit of very strong Coulomb repulsion 
the spin degrees of freedom are commonly modelled by 
the spin-1/2 Heisenberg chain. This has proved to pro- 
vide an adequate description of recent inelastic neutron 
scattering experiments on SrCuC-2 00- In this limit the 
low energy dynamics is dominated by spin flips generated 
through virtual hopping processes of electrons. However, 
for general Mott insulating materials there is no reason 
for the Coulomb repulsion to be much stronger than the 
electron hopping amplitude. Indeed, this situation is re- 
alized in the Bechgaard salts. A natural question to ask 
then is how electron itineracy affects the spin dynamics. 
This issue has recently attracted much attention in the 
context of the two-dimensional Mott insulating cuprate 
La2Cu04 0]. The experimentally observed spin wave 
dispersion was found to depart significantly from the 
Heisenberg form, revealing the presence of ring-exchange 
interactions. This is a direct manifestation of electron 
itineracy in the sense that ring-exchange interactions are 
generated by higher order virtual electron hopping pro- 
cesses. 



Motivated by these findings we investigate the effects of 
electron itineracy on the unusual spin correlations in Q1D 
Mott insulators. We do so in the simplest model that 
incorporates the relevant physics, the one dimensional 
half filled Hubbard model. The Hamiltonian is given by 



t, a 

+u^i,r-b(n l ,i-b, (1) 



where cf a creates an electron with spin a =f, I at site I, 
n^o- = c^Q cr, and U is the on-site Coulomb repulsion. In 
momentum space, the kinetic term yields a cosine band, 
e{k) = — 2icos(fca ), where a is the lattice spacing. At 
half- filling, kp = ±h/2ciq. The elementary excitations of 
this model are termed spinons and holons. The spinons 
are gapless, uncharged, and carry spin 1/2. Their energy 
and momentum are conveniently parameterized in terms 
of the spectral parameter A 0] : 
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where Jq(x) and J\{x) are Bessel functions, and E s (0) 
plays the role of a "spinon bandwidth" . The holon and 
antiholon, are gapped, spinless, and carry charge |e| and 
— |e| respectively. The "charge gap" A, defined as the 
minimum of the holon energy, increases with U and is 
given by: 
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These excitations form a basis of scattering states in 
which one may explore the dynamical response of ID 
Mott insulators. 

As usual, the dynamical structure factor is obtained 
from the imaginary part of the Fourier transform of the 
spin-spin correlation function. In view of the SU(2) spin 
symmetry of the Hamiltonian we consider: 



S zz (Lu,q) = -Im (xjsfw.ffJIo-Hj-iw), 



(5) 
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where r\ — > + and the susceptibility (in Euclidean time) 
is given by 



dx dTe^-^S'^T^S'^O)). (6) 



The spin operator is given by S z = (n-f — njJ/2. In the 
first instance, it is instructive to study the known limits 
of large and small U/t. In the limit U ^> t, the half 
filled Hubbard model reduces to an isotropic Heiscnbcrg 
antiferromagnet with exchange J ~ 4t 2 /U In this 
case, the structure factor is well approximated by the 
so-called Muller ansatz 
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Whilst more detailed results can been obtained using in- 
tegrability |lOj ] this approximation is sufficient for our 
current discussion. The lower and upper spinon bound- 
aries, ujl = (tt/2) J I singao| and u>u = ttJ | sin(gao/2)|, 
follow from the spinon dispersion relations @ and 
These boundaries delimit a continuum of two spinon ex- 
citations, and we plot this in Fig. |l(a)| It is readily seen 
that the spectral weight (plotted darker) is concentrated 
on the lower spinon boundary. The power law singular- 
ity extends down to lu = in the vicinity of q = tt, and 
reflects the gapless nature of spinons. 

In the limit of large external freqencies and U <C t, one 
may treat the Hubbard interaction perturbatively. Using 
the tightbinding propagator Q(u>,q) = clq [iu> + 2t cos q] _1 , 
a summation over "bubble" diagrams yields: 
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where the U = single bubble is given by 
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(9) 
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We plot the corresponding structure factor in Fig 
In this limit, significant spectral weight is concentrated 
on the upper spinon boundary. (It is evident that this 
large to approximation fails to capture the low frequency 
spinon divergence; see discussion below.) Thus apart 
from a rescaling of the spinon bandwidth, to which we 
shall return below and in Fig. there is a manifest redis- 
tribution of spectral weight as one changes the ratio U/t. 
It is therefore highly desirable to build up a better pic- 
ture for what happens in the intermediate regime U ~ t. 
We address this pertinent issue using the techniques of 
integrable field theory. 

At low energies, we may linearize the non-interacting 
fermionic spectrum around ±|fc.p|. The resulting model 
may be bosonized in terms of two Bose fields, $ s and 
<f> c . These correspond to the collective spin and charge 
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FIG. 1: Density plots of the dynamical structure factor: (a) 
U 3> t Muller approximation, (6) U <C t result for U/t = 
0.4. At large (small) U/t the spectral weight is on the lower 
(upper) spinon boundary. 



degrees of freedom respectively. The continuum Hamil- 
tonian density separates in two commuting pieces: 



where 



H = H s +H C 



g : J J :, 
g : I I : . 



Here J and I are the SU(2) spin and isospin currents, 
v s = VF — Uao/2ir, v c = vp + Uao/2w and g = 2Uclq. The 
interaction in the charge sector is marginally relevant and 
leads to the formation of the charge gap. The interaction 
in the spin sector is marginally irrelevant and will be 
neglected. This description is valid for A <C 4t and for 
U ;$ 3i. In order to enhance the clarity of presentation, 
we will mainly work in the so-called scaling limit: 
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U/t -> 0, A 



me -2nt/U fixe(L 



In this limit v s coincides with v c and Lorentz invariancc 
emerges. The charge sector is governed by the SU(2) 
invariant sine Gordon model, and the spin sector by a 
massless Gaussian model. This simplifies the ensuing 
formulae, whilst capturing the relevant physics. Where 
appropriate, we give the results for v s =/= v c . 

In this low energy description, the z component of the 
spin operator takes the form: 

S* = -^a 8 $ 8 + (-l)*^cos^sm%. (10) 

47T TT I I 

We draw attention to the fact that both the "spin" and 
the "charge" fields enter the staggered component of the 
spin operator. Excitations of the charge field, and thus 
electron itineracy, may therefore influence magnetic cor- 
relations in the vicinity of q ~ tt. In the first, crude 
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approximation, the operator cos is replaced by its 
constant vacuum expectation value ( | cos ^f - | ) , and 
correlators are given by the critical spin sector alone. 
However, this is merely the leading term in a system- 
atic expansion of the charge sector. In contrast to naive 
expectations, these terms are shown to be significant. In 
the subsequent analysis, we shall calculate the "itiner- 
acy corrections" to the magnetic structure factor arising 
from these charge excitations. We employ a number of 
results developed in the study of density-density correla- 
tions 0|- A s follows from equation (|10|l the staggered 
part of the spin-spin correlation function is given by 

{S z (r,x)S z (0,0)) stagg . = (-1)* ^M, (11) 

2tt z vf r 



where r = a/t 2 + x 2 jvp and we have defined the contri- 
bution from the gapped charge sector 

^(r) = (0|cos^(r,a:)coBy(0,0)|0). (12) 

We introduce a basis of states | 9 n . . . 9\ ) e „...ei where the 
index e = ± denotes holon and antiholon, and 9 are their 
rapidities. Inserting the resolution of the identity 

1 = |0)(0| + 

EE /" ••*>«....« —Ho n ...e 1 \ 



n—l £i 
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one may develop i|12|> in a systematic expansion over 
multi-particle intermediate states: 



T(r) = T (r) + Mr) + 



(13) 



The leading term is obtained from the vacuum expecta- 
tion value of the operator 
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S -|C| S 
v F 



(14) 



and is a constant. The dimension of (0|cos^|0) is 



[L]^ 1 / 2 . The dimensionless constant C is pinne d by the 
short distance properties of the charge sector |l3j and 
fixes the normalization of the structure factor. That is 
to say, to leading order, the staggered spin-spin correla- 
tion function is approximated by the conformal 1/r de- 
cay. Since the operator cos $ c /2 does not couple to single 
holon states the next contribution may be written: 
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x exp [— m(c\\9i + c1i#2)|t| — im(s\\9i + s\i92)x/vf] ■ 
Defining 9 = 9\ — 9 2 , the form factor is given by [Til . [l5| 



f cos- {0) + _ = (o|cos-^|0 3 0i>+_ 



(15) 



E(9) = exp 



dx sin 2 ([# + in]x/ir)e 



i x sh(2x)ch(a;) / 
In view of (|13f) we may expand the structure factor 



S**(u,q) = ^TSI*(u;,q) , 



(16) 



(17) 



where S 2 *(uj 7 q) denotes the contribution due to inter- 
mediate states with 2n (anti)holons. The "conformal" 
contribution is easily evaluated: 
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where Q(y) = 1 for y > is the Heaviside step function, 
and s = \J uj 2 — vp q 2 . Performing the sum over isotopic 
indices and the necessary integrals the "itineracy 
correction" is found to be 

S?{u,± + = 

ao vf 27r z 



d6 8 sh^ 
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where 9' = arch (s/2m), 7 = 2ch9. We plot these con- 
tributions to the structure factor in Fig. |2] For v s 7^ v c 
(|19|l may be generalized to 
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d9d& |/< 



(29), 



e(n+)e(n_), 



where f2± = ui±v s q-2A ch9 (ch.9'±(v s /v c ) sh.9'). In the 
vicinity of the threshold (s = 2m), the integral in l|19(l 
may be obtained by Taylor expanding the integrand. The 
leading threshold behavior as s — > 2m + is 
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(20) 



In contrast to the conformal contribution, this term goes 
to zero at its threshold and is free of singularities. At 
high energies s 2m, S 2 z (uj, +q) — > constant. In real 
space, the asymptotics of this correction (evaluated in 
the saddle point approximation) yield: 
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for (mr) 3> 1, where 
2V2 



A 



■ exp 



dx e x sh x 
x ch 2 x 



(21) 



(22) 



It is apparent from Fig. [21 that the itineracy cor- 
rection becomes significant for s ~ 10m; for v s =/= v c , 
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FIG. 2: The staggered dynamical structure factor. We plot 
the conformal contribution, the holon-antiholon "itineracy 
correction" and their sum. The inset shows the ratio of the 
charge gap to the spinon bandwidth as a function of U/t. 
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FIG. 3: Approximate collapse of the rescaled spinon con- 
tinuum. We plot the rescaled spinon energy E' a (A) = 
E a (A) / E a (0) versus momentum P a (A). The arrows show the 
transfer of spectral weight as U/t is increased. 



s ~ 10A. In order to observe this feature within the 
two spinon continuum, we seek to maximize A subject 
to 10 A < 2E S (Q) < 40t. As may be seen from the inset 
of Fig. [21 this condition is satisfied for 2 < U/t < 3. In 
particular, the spinon bandwidth E s (0) decreases with 
increasing U/t. In view of this we rescale the spinon en- 
ergies by this factor. As shown in Fig. |3 the rescaled 
dispersion relations collapse (approximately, but remark- 
ably well) onto a single curve. We may therefore compare 
systems with different values of U/t on a single plot. In 
conjunction with the limiting cases, [/ « t and U 3> t, 



our detailed analytic results support the simple yet com- 
pelling view depicted in Fig. |3| As U/t increases spectral 
weight is transferred from the upper rescaled boundary, 
and builds up in both the conformal divergence and the 
itineracy correction. The downward shift at high fre- 
quencies is evident from the bubble summation |JSJ and 
the build up at low frequencies originates in the increase 
of the charge gap. 

In this letter we have determined dynamical spin cor- 
relations in ID Mott insulators. Away from the well- 
understood "Heisenberg limit" of very large Mott gaps, 
electron itineracy effects are shown to be quite significant. 
On the basis of our calculations we have put forward 
a simple picture describing the spectral weight transfer 
in the dynamical structure factor as a function of the 
strength of the Coulomb repulsion. We expect our find- 
ings to be of relevance for inelastic neutron scattering 
experiments on quasi-lD small-gap Mott insulators. 
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